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ABSTRACT 

We examine the global stability of an exponential stellar disk embedded in a dark mat- 
ter halo and constrain the mass-to- light ratio of observed disks in the I-band, T/. As- 
suming only that the radial surface density distribution of disks is exponential, we de- 
rive an analytic upper limit: T/ ^ 2.4/i, for a Hubble constant of lOO/i km s~"'^Mpc^^. 
Using iV-body simulations we derive a stability criterion significantly different from 
that of previous authors. Using this criterion we argue that, almost independent of 
the concentration of the halo, T/ ^ 1.9h. We discuss this result in relation to other 
independent determinations of T/, its limitations and its implications for theories of 
disk formation and barred galaxies. 
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1 INTRODUCTION 

Stellar disks are a potentially powerful probe of the mass 
distribution in spiral galaxies. In particular, the rotation 
curve (defined in Section |^ is a direct probe of the grav- 
itational force in the disk. The fact that the rotation curves 
of spiral galaxies are rather flat is usually taken to imply 
the presence of an extended halo of dark matter (e.g. Free- 
man 1970, Persic & Salucci 1991). Dark halos have been 
traditionally modelled by isothermal spheres with homoge- 
neous cores, and assigned a mass-to-light ratio according to 
the 'maximum-disk' hypothesis, in which the largest mass 
possible is assigned to the disk consistent with the rota- 
tion curve (Carignan & Freeman 1985). Recent theoretical 
work has suggested that the isothermal form of the halo may 
not be appropriate (Navarro, Frenk & White 1996, Moore 
etal. 1997, Kravstov etal. 1997). Rix & Courteau (1997) 
have recently argued that the TuUy-Fisher relation (TuUy & 
Fisher 1977) has too little scatter to be consistent with the 
maximum-disk hypothesis. 

Dalcanton, Spergel & Summers (1997), and Mo, Mao & 
White (1997) point up the need for a consistent and accurate 
picture of disk mass-to-light ratios for the theory of disk for- 
mation, and to understand the TuUy-Fisher relation. As was 
pointed out by Rix & Courteau (1997), the maximum disk 
hypothesis and the TuUy-Fisher relation are together incom- 
patible with a near universal mass-to-light ratio, T. Yet, the 
stellar populations of spiral disks at z = are generally be- 
lieved to have rather uniform properties. In particular they 
all have similar colours, which leads one to believe that they 
should have similar T. 
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Ostriker & Peebles (1973) argued that extended halos 
of dark matter in disk galaxies are required to stabilize the 
disks against global bar instabilities. They were able to draw 
very strong conclusions from A-body simulations with only 
300 particles. Efstathiou, Lake & Negroponte (1982) (ELN) 
conducted larger A'^-body experiments specifically with disk 
components appropriate for real galaxies, and concluded 
that Tboi < 1.5 ±0.2. 

ELN compared their numerical results with only 12 
galaxies due to limited data at that time. Since then, ob- 
servations of spiral galaxies have improved tremendously 
owing to the great interest in the TuUy-Fisher relation as 
a distance indicator (Giovanelli etal. 1997 and references 
therein). Large samples of spiral galaxies are now available 
(e.g. Mathewson & Ford 1996; Courteau 1996, 1997). In this 
paper, we re-examine the constraints from considerations of 
global stability which may be placed on the mass-to-light 
ratio of disks. 

In the next section we review the properties of exponen- 
tial disks, and their relation to real disk galaxies. In Section 
^ we review the physics of global instability in disks. In Sec- 
tion^ we describe the observations of disk galaxies, and how 
they can be used to constrain T. In Section H we describe our 
A-body simulations. In Section ^ we discuss the results of 
independent determinations of mass-to-light ratios of galac- 
tic disks. In Section ^ we discuss the implications of our 
results and draw conclusions. 



2 EXPONENTIAL DISKS 

The luminous disks of spiral galaxies are commonly mod- 
elled by an exponential surface brightness distribution: 
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where R is the usual cyhndrical radius, and La is the total 
luminosity of the disk. Here we collect some notation and a 
number of useful results relating to exponential disks. 

The disk has a mass Md, and a mass-to-light ratio in 
the /-band T/ in solar units. Thus the surface mass density 
of the disk is 



(2) 



The gravitational potential in the disk $ is conveniently 
decomposed into contributions from the disk and a halo: 



$ = $d + <E>h. 



(3) 



(We use the subscripts 'd' for 'disk', and 'h' for 'halo' 
throughout.) We assume for the present purposes that the 
halo is spherical, and usually we think of it as being com- 
posed of dark matter, but it may also contain a stellar com- 
ponent (e.g. the 'bulge' of an earlier type spiral). 

The speed of test particles on circular orbits Vc as a 
function of R is given by 

^UR)--^. (4) 

For the sake of definiteness, we shall refer to Vc{R) as the 
'true rotation curve' of the system (or just 'the rotation 
curve'). The observed line-of-sight velocity of a tracer pop- 
ulation corrected for inclination we refer to as the rota- 
tion curve of the tracer. For instance, the HI rotation curve 
(apart from small contributions from turbulent motion) is 
thought to be a good measure of the true rotation curve as 
long as the system is axisymmetric. Some of the samples of 
galaxies used in this paper have rotation curves measured in 
Ha, but analyzed in a way which is supposed to maximize 
the agreement with HI (Raychaudhury et al 1997, Courteau 
1997). 

Following Efstathiou, Lake & Negroponte (1982) we de- 
fine the dimensionless quantity 

(5) 



(GMd/i?d)i/2' 

where Vm is the maximum value of Vc, and G is the gravita- 
tional constant. The rotation curve of an isolated exponen- 
tial disk ($h = 0) is (Freeman 1970): 



y\h{y)Ka{y) - h{,y)K^{.y)].. 



(6) 



2 2GMd 2 

'-'d(^) = —5- " 

where y = R/{2Rd). h and Ki are modified Bessel functions. 
An isolated disk has em = fd ~ 0.63 and a disk embedded 
in a halo has em > fd- 



exemplified by the Maclaurin sequence of gaseous spheroids 
(Chandrasekhar 1969): near spherical members of the fam- 
ily are stable; they become first secularly and then dynami- 
cally unstable as the eccentricity is increased (see for exam- 
ple Christodoulou, Schlosman & Tohline 1995). A similar 
pattern is found for the stellar analogues of the Maclaurin 
spheroids. Stability may also be conferred by adding an ex- 
ternal potential — in the extreme limit that the disk is com- 
posed of test particles it is trivially stable. As was pointed 
out by Ostriker & Peebles (1973) this mechanism can be 
used on rather general grounds to argue that spiral galaxies 
contain halos of dark matter. 

There has been some debate as whether a disk with a 
given eccentricity and/or external potential is expected to 
be stable on theoretical grounds. Many attempts have been 
made to find a single parameter which delineates a boundary 
between stability and instability. Ostriker & Peebles (1973) 
used as a parameter the ratio 



top = rmoan/W 



(7) 



of the mean kinetic to total potential energy. The stellar 
disks studied by Ostriker & Peebles (1973) were stable if 
top < 0.14 to very high accuracy. Remarkably, this criterion 
was shown to unify their results with all existing results on 
the stability of disks. Subsequent authors pointed to various 
insufficiencies of top as applied to specific cases (Miller 1978, 
Miller & Smith 1979). Efstathiou, Lake & Negroponte (1982) 
refined the definition to patch up some of these difficulties, 
defining a related parameter t*, but in the end concluded 
that em (cf. equation |^ was a better indicator of stability 
for their models. In particular they found that all models 
with em ^ 1.1 were unstable, whereas the stability boundary 
in t* depended on the halo concentration. Christodoulou, 
Schlosman & Tohline (1995) proposed a further refinement 
of t* by analogy with the Maclaurin series, defining a new 
parameter 

a = yi/i' (8) 

where t = T /W and / is a form factor dependent on the 
shape of the disk. Both t and / were given precise definitions 
only for the Maclaurin sequences. They showed that a led to 
a stability criterion which was broadly consistent with that 
of ELN. However, a was not defined properly for disks with 
external potentials, and they failed to address the issue of 
the halo concentration. 

The importance of global stability is this: the disks in 
spiral galaxies appear to be long lived, and the majority do 
not contain bars. In the next section we describe how the 
apparent stability of galactic disks can be used to constrain 
their mass-to-light ratio. 



3 GLOBAL STABILITY 

An isolated thin exponential disk is known to be violently 
unstable, as are all isolated thin disks. The fastest grow- 
ing instability is generically global (the whole system is 
deformed) and bar shaped. The instability in conservative 
and/or coUisionless systems is generally purely dynamical 
(growth timescale of order the crossing time). Secular in- 
stabilities set in in systems with dissipation and/or colli- 
sions generally on longer timescales. The disk may be sta- 
bilized if the disk is thickened in the vertical direction, as 



4 COMPARISON WITH OBSERVATIONS 

To compare a theoretical instability criterion with observa- 
tions, we need both detailed rotation curves (to obtain Um) 
and surface photometry (to obtain Rd)- We have examined 
the data from a wide variety of sources. The largest data 
set is that of Mathewson & Ford (1996) (MF) which has 
rotation velocities and /-band photometry for a sample of 
nearly 2500 Southern spiral galaxies taken from the ESO- 
Uppsala catalogue. The majority are relatively late types: 
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Figure 1. The quantity eih~^/^ (as defined in equation is 
shown versus 7-band central surface brightness for the sample 
of spiral galaxies of Mathewson & Ford (1996) (with crude bulge 
correction, see text). On this plot, a fixed value of tm corresponds 
to a horizontal line with amplitude e-m^'^ ih~^)^/^ . The horizontal 
solid line shows the value of em for a self-gravitating disk (em = 
0.63) with T I = 2.14/1. The two shaded lines bracket the predicted 
ranges (see text). Barred galaxies are shown as triangles. 



Table 1. Observational data and summary of limits on Tj. Col- 
umn (1), Name of data set; (2), number of galaxies; (3) median 
central surface brightness (/xo) in mag/sec^; (4), e; minimum, 
maximum, median and 10% quantile values; (5) Tj limit de- 
rived from 10% quantile and assuming em ^ 0.63; (6) assuming 
em > 0.72. See text for description of data set names. The row la- 
belled 'Air contains the combined properties of all the preceding 
rows. MFB refers to the subset of MF which is barred. 



Set 


N 


(w> 


e; [min,max,med,10%] 


Ti(l) 


Ti(2) 


MF 


2446 


19.97 


[0.43,7.50,1.19,0.92] 


2.14 


1.64 


C97 


304 


20.14 


[0.54,2.03,1.13,0.90] 


2.02 


1.55 


C96 


316 


20.39 


[0.76,1.84,1.15,0.96] 


2.34 


1.79 


SML 


23 


22.87 


[0.54,5.72, 1.95, 1.04] 


2.72 


2.08 


SMH 


23 


20.97 


[0.67,3.56,1.32,0.80] 


1.62 


1.24 


BBFN 


371 


18.47 


[0.35,2.10,0.99,0.70] 


1.24 


0.95 


DG 


84 


20.05 


[0.59,1.77,1.04,0.84] 


1.79 


1.37 


RBKG 


25 


20.89 


[0.89,1.78,1.28,1.08] 


2.93 


2.24 


B 


12 


19.89 


[0.67,1.96,1.12,0.84] 


1.77 


1.35 


All 


3604 


20.01 


[0.35,7.50,1.16,0.88] 


1.96 


1.50 


MFD 


2446 


19.97 


[0.45,7.50,1.30,0.98] 


2.41 


1.85 


MFB 


175 


19.94 


[0.54,4.37,1.20,0.90] 


2.05 


1.57 



of those 2275 for which Hubble types are given, 1055 are 
Sc or Sbc; 814 are Sb; and 5 are Sa. We convert the pub- 
lished photometric quantities to and the central surface 
brightness by assuming an exponential profile. Details 
are given in an Appendix. The Appendix also describes how 
we make a crude bulge subtraction from the MF data, the 
results being labelled MFD in Table ^. Since Mathewson & 
Ford (1996) provide the largest uniform data set we adopt 
it whenever we quote a numerical result (including bulge 
subtraction when not explicitly specified). We list the other 
data sets for comparison purposes. 



The data of Courteau (1996) and Courteau (1997) are 
in principle the same, but one paper gives kinematic mea- 
surements and the other photometric, and we list them sep- 
arately. To obtain the missing data in these cases we use 
the observed TuUy-Fisher relation (Giovanelli et al. 1997; 
Shanks 1997): 

A1/-51og/i = -(21.00±0.02)-(7.68±0.13)(log W^-2.5),(9) 

where h is the Hubble constant in units of 100 km s~^Mpc, 
and where W is the inclination-corrected width of the HI 
line profile. The maximum of the rotation curve Um is given 
by I'm = W/2. Courteau (1997) has measured rotation not 
in HI, but in H^. Raychaudhury etal. (1997) show that to 
a good approximation values of W determined from Hq are 
consistent with those determined from HI. 

For data sets without /-band quantities we convert 
using for all galaxies B — I = 1.7 (de Jong 1996) and 
r - I = 0.77 (Rix & Courteau 1997). When or /xo are 
missing we calculate them assuming an exponential profile. 

The other data sources are as follows with a short de- 
scription of the quantities provided. C97: (i?d, fm) provided 
by the author from Courteau (1997). C96: Courteau (1996) 
{Rd, r-magnitudes). SML: de Blok & McGaugh (1996) low 
surface brightness galaxies (i?d, «m, B- magnitudes). SMH: 
de Blok & McGaugh (1996) high surface brightness galaxies 
(ditto). BBFN: Burstein etal (1996) spiral types {vm, effec- 
tive B-band photometric quantities). DG: Dale etal (1995), 
(i?d, Wm, /-magnitudes and ^o)- RBKG: Raychaudhury etal 
(1997), {Rd, Wm, /-magnitudes). B: Bottema (1993), {Rd, 
Wm, iJ-magnitudes) . 

The luminosity (in solar units) is derived from listed 
magnitudes using Lj = iq0-4(415-a1/) then directly 

compute the quantity 

(10) 



{GLi/RdY'^' 
which is related to tm by 

ei = euilj ■ 



(11) 



The quantity ef has the units of a mass-to-light ratio, and 
indeed it is a measure of the total mass (including dark mat- 
ter halo) contributing to the rotation curve. Let us define a 
quantity 



Tf\R) = 



Vc{R)R 
GLi{R)' 



(12) 



which measures the total mass-to-light ratio as a function 
of radius. For an isolated disk it is a constant (= T/) and 
with an extended dark halo it increases with radius. The 
maximum rotation velocity in general occurs &t R = /Zm > 
Rd, and the luminosity enclosed is L/_m < Li, hence ef < 
Tf*(/?m). For an isolated disk /Zm = 2.2/Zd and L/,m = 
0.65L/, so T/ = T*/*(/Zm) = 3.4ef . 

In Figure |l], we show ei as a function of central sur- 
face density for the MFD data. The extra factor of h 
in the abscissa makes the plotted quantities independent of 
the Hubble constant. The figure reveals a marked correla- 
tion between e; and ^o- This is expected if T/ is independent 
of ^0 in the standard picture of disk formation (Fall & Ef- 
stathiou 1980). The trend can be understood with a simple 
model where the dark matter halo is assumed to be a singu- 
lar isothermal sphere and disk self-gravity is ignored (Mo, 
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Mao & White 1997). In this model, the disk scale length and 
surface mass density scale as 

7?dOcAiih, Eo oc ^ oc (13) 

where is the fraction of halo mass that settles into the 
disk, A is the dimensionless spin parameter, Vm is the halo 
circular velocity and H{z) is the Hubble constant at redshift 
z (see Mo, Mao & White 1997 for details). The TuUy-Fisher 
relation is given by 

T A 3 A _ _r Ho , . 

Lj=Av,^, A = /,^-^^, (14) 

where /; is a dimensionless constant which depends on the 
cosmology and halo profile (Mo, Mao & White 1997). From 
equations ( p^ & have 

7,-1/2 ^1/2 .-1/2 f Ho \ 

e^h / ocA/ A (^-^j , (15) 
and 

where ^o = So/T/. Notice that m^, Tj and the Hubble 
constant come into the above expressions only through A 
defined in equation (|l^). From equations ( |l^ ) and ( [iq ) we 
obtain the dependence of e; on as 

eih-'^' oc f,-'^\K'A-'^\ (17) 

For a given central surface brightness, the scatter in e; is 
determined by the range in circular velocities and the scat- 
ter in the TuUy-Fisher amplitude A. Since the scatter in 
A is approximately a factor of 2 and the velocity range 
is ~ 100 — 300 km s~^, we expect a scatter of roughly 60 
percent. In Figure 0, we overlay the predicted slope and 
scatter on top of the data points with the normalization 
chosen to reproduce the observed median value of e; in the 
range fio £ (20,21). The lines are derived from the 25 and 
75% quantiles of A/v„^ of the data — the equivalent scat- 
ter in A is a factor of ~ 1.9 or 0.7 magnitudes. As can be 
seen, the predicted slope and scatter are consistent with the 
observed data points. Thus, the data are consistent with 
the assumption that the value of A is independent of /io 
for these galaxies. At the fainter end the observed e; is 
slightly higher than the model prediction, implying either 
Tj is higher or md is lower for low-surface-brightness galax- 
ies. The former would result from a lower star formation 
efficiency, as suggested by observations (e.g. McGaugh & 
de Blok 1997). Since e; oc A^/^ (cf. equation |l|), a sub- 
stantial angular momentum loss would lower e; significantly 
and make many disks unstable. Hence, the stability of disks 
require approximate conservation of angular momentum in 
the disk formation process. The normalization in equation 
( p^ is dependent on the detailed halo profiles and other 
factors. A detailed comparison with the observations would 
also need to take into account observational biases (e.g, in 
surface brightness), which is beyond the scope of the present 
work. 

Assuming that the disks are exponential already sets a 
limit on T/: they should all have Em > 0.63, the value for an 
isolated disk. To be conservative, let us suppose that 90% 



of galaxies should have em > 0.63, and ascribe the remain- 
ing 10% to deviations from exponential and/or measurement 
errors. In Figure ^ the solid horizontal line marks the 10% 
quantile of ei in the data of Mathewson & Ford (1996) (after 
crude bulge subtraction, see Appendix). This corresponds to 
em — 0.63 for T/ = 2.41/i. When the bulge subtraction is 
not carried out the implied value of T/ goes down by about 
20 percent to T/ — 2.14/i. The quantiles and implied limits 
on Ti for the other data sets are given in Table 0. The dif- 
ferences between the data sets are on the whole not large. 
Burstein etal (1996) (BBFN) give only effective quantities, 
and hence we believe this data is more heavily contami- 
nated with bulge contributions than the others. This would 
account for the lower values of e; and T/ (and also for the 
brighter median /io). 

If Tj > 2.41/i more than 10% of galaxies would have 
em < 0.63 which would be incompatible with their hav- 
ing exponential surface density distributions. Thus a con- 
servative upper limit on the mass-to-light ratio of disks is 
T/ < 2Alh. Note that this limit applies only under the as- 
sumption that Tj is independent of fio- This assumption is, 
as we argued, consistent with the observational data, if we 
adopt the standard picture of disk formation. The derived 
limit is also consistent with independent measurements of 
T/ (see Section 6 for details). However, since e; is higher 
for low-surface-brightness galaxies, a higher upper limit on 
Tj is still allowed for these galaxies without violating the 
constraint em > 0.63. 

A tighter upper limit on T/ is provided by demand- 
ing that disks are dynamically stable. The stability criterion 
proposed by ELN (em < 1.1) implies a very low value of T/: 
the 10% quantile giving an upper limit of T/ < 0.79h for 
the MFD data. This upper limit is very low and is incon- 
sistent with direct measurements of Tj. To try and resolve 
this conflict we carried out our own A'^-body simulations to 
locate the stability boundary independently. 



5 NUMERICAL SIMULATIONS 

Here we describe A'^-body simulations of exponential disks 
embedded in dark halos. The halos we use, in common with 
ELN, are rigid (i.e. they are modelled by a fixed poten- 
tial). The halo density profile used was either a Hernquist 
(1990) profile, or a truncated isothermal sphere as described 
in Hernquist (1993). A comparison between the two halo 
types probes the effect of halo concentration on disk stabil- 
ity. The disk is locally isothermal: 

p{R,z)^^^sech^iz/zo) (18) 

ZZo 

with a constant vertical scale height of zo — 0.2i?d at all 
radii. The radial dispersion velocity is given by 

(4) = Cexp(-i?/i?d) (19) 

with C constant. This is consistent with observations of 
galactic disks (van der Kruit & Searle 1981, Bottema 1987). 
The constant C is chosen by fixing Q — 1.2 (Toomre 1963) 
at R — 2ARd (this is then roughly the minimum Q in the 
disk). Galactic stellar disks may well have larger values of 
Q (Lewis & Freeman 1989), but this would lead to greater 
stability, and hence the limit we derive will be conservative. 
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Table 2. Simulation parameters. Lengths in units of R^; masses 
in units of M^. LH=Hernquist profile; IS=isothermal. Afh is total 
halo mass. The IS halos were truncated exponentially at r = 8 
(TIernquist 1993); a is halo scale length; is defined in equation 



Name 


Halo 


Mh 


a 


Em 


bar 


1001 


LH 


1 


1 


0.78 


no 


1002 


LH 


2 


1 


0.92 


no 


1004 


LH 


4 


1 


1.15 


no 


1005 


LH 


.5 


1 


0.70 


yes 


1101 


IS 


2.5 


1 


0.78 


no 


1102 


IS 


5 


1 


0.92 


no 


1104 


IS 


10 


1 


1.16 


no 


1103 


IS 


2 


1 


0.75 


no 


1105 


IS 


1.25 


1 


0.70 


yes 




Figure 3. The endpoints of two of the simulations (showing only 
4096 particles to avoid crowding). Simulation 1005 clearly shows 
a strong bar, while 1001 does not. The rotation curves of these 
two models are shown in Figure ^. The figures are bounded by 
the simulation box which is SR^ on a side. 




10.0 - 
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Figure 2. The rotation curves of some of the models in Table pj^^^^ ^ ^j^^ strength of the bar |Ao| as a function of time 

|5|. The lower curve is that of an isolated exponential disk. simulations in Table | Simulations 1005 and 1105 show 

clearly the presence of a bar, as shown in Figure ^. The time unit 

is Td = 27r(R3/GMd)l/2. 



The disk equilibria were set up as described by Hernquist 
(1993) using a code provided by Professor Hernquist himself. 

We used the AP^M code of Couchman, Pearce and 
Thomas (1995) compiled for isolated (not cosmological) ini- 
tial conditions, and modified to include the rigid halo con- 
tribution. The simulations had 16384 disk particles, and a 
softening length of 0.05i?d. Each simulation was run for ap- 
proximately 7 times l-K^Rl/GMd)^^^ ■ Table | lists the sim- 
ulation parameters in full. The rotation curves of some of 
the models are shown in Figure |^. The Hernquist halos have 
cuspy profiles and hence the rotation curves rise much more 
steeply than in the isothermal case for the same scale length. 
The Hernquist rotation curves are also more peaked than the 
isothermal case. Peakiness is not in this instance an indica- 
tion of low em but rather of the finite mass of the halo: the 
rotation curves are roughly Keplerian for 3. 

We can evaluate the strength of the instability visually 
and quantitatively. The endpoints of simulations 1005 and 
1001 are shown in Figure ^, and a bar is clearly visible in 
the former. As an objective measure we use a logarithmic 
spiral decomposition of the surface density (Anathassoula & 
Sellwood 1986). We define the quantity 



^p-^Yl ""^ exp[i(2ej +plnRj)], (20) 

where p is the logarithmic spiral wavenumber, irij is the 
mass, and {Rj,9j) are the cylindrical co-ordinates of a par- 
ticle labelled by an index j. The factor outside the summa- 
tion scales Ap so that the expected value is unity if there is 
no significant structure at a given p. The factor of 2 in the 
argument of the exponential means that we pick out density 
perturbations with m = 2 symmetry. For p > we have a 
trailing spiral and for p < a leading spiral. A bar has p — 0. 
We measure Ap approximately 40 times during each simula- 
tion, and consider the quantity |j4o| defined as the average 
over the range p G [—1,1]. In Figure ^ we show \Ao\ as a 
function of time for all the simulations. Those with visually 
identifiable bars (1005, 1105) show up clearly as having a 
significantly larger value of \Ao\. Thus both methods agree. 
We give an indication of whether a bar was present or not 
in Table |. 

It follows that the boundary ec between global stability 
and a bar-forming instability lies between em — 0.70 and 
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Em = 0.75. The 10% quantile in Figure |l] combined with 
these values of tc gives an upper limit to the mass-to-light 
ratio Ti < (1.85±.l)/i. The formal error only reflects the un- 
certainty in the value of ec . The upper limit would be higher 
if one had more tolerance for disks with apparent tm < Ec. 
Nevertheless the limit is in agreement with independent di- 
rect measurements (see Section ^ . 

The critical value ec is much lower than that obtained 
by ELN for their 'FE' models (ec « 1.1). ELN also con- 
ducted simulations with isothermal halos similar to ours, 
but they did not cover as much parameter space as their FE 
models, and were not analyzed in much detail. The obvious 
differences between our simulations and those of ELN are 
that our simulations are 3-dimensional, and that they have 
slightly better formal resolution. Resolution would not seem 
to explain our different conclusion. Indeed, we carried out 
simulations with A'' = 1024 and a significantly larger soften- 
ing length and reached the same conclusion. The difference 
could arise from the flnite thickness of our disks. They have 
a vertical scale height of 0.2_Rd, implying a formal eccen- 
tricity of e = y/1 — .2^ = 0.97. Isolated stellar Maclaurin 
spheroids are already unstable at e = 0.97. The added ex- 
ternal potential may be enough to stabilize our disks but 
not those of ELN with e = 1. A further difference between 
our simulations and those of ELN is the value of Toomre's 
Q and its dependence on radius. Most of the simulations of 
ELN use Q{R) — constant = 1.05, whereas we use a profile 
with a broad minimum of Q = 1.2 at ii ~ 2.4iid. ELN also 
ran a few simulations with Q a decreasing function of radius, 
and reported that their results were unchanged. The lower 
value of Q would make their disks more prone to bar for- 
mation. Note however that typical values for galactic disks 
are even higher than in our simulations (Lewis & Freeman 
1989, Bottema 1987, von Linden & Fuchs 1997). 

Our results support the conclusion of ELN that is 
a good indication of stability, independent of halo type and 
concentration (the Hernquist halos are more concentrated 
than the isothermal ones). One implication is that a modi- 
fied Ostriker & Peebles criterion such as that considered by 
Christodoulou, Schlosman & Tohline (1995) is unlikely to 
be universal. 



6 INDEPENDENT DETERMINATIONS OF T 

Limits on the mass-to-light ratio of the Galactic disk in 
the solar neighbourhood can be derived from a combina- 
tion of kinematic measurements and star counts. Kuijken 
& Gilmore (1989) derive a local surface mass density in the 
disk of 40Mq /pc^ , and star counts give a 1/-band luminosity 
density of 15Lq/pc^ (Gould, BahcaU & Flynn 1996). Divid- 
ing mass by light we obtain Ty = 2.67, and thus T/ ~ 1.9 
(assuming V — I = 1.0). This number is independent of 
h, and hence is rather high compared with the upper limit 
derived from the A^-body simulations. The data quoted by 
Bottema (1993) give a value of ei = 1.37 for the Galaxy, 
so for Ti « 1.9 we obtain em ~ 1.0. Thus despite having 
high Ti we conclude that the disk of the Galaxy should be 
globally stable. 

Mass-to-light ratios can also be derived from pure stel- 
lar population synthesis arguments, although there is al- 
ways some uncertainty arising from the poorly known initial 
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Figure 5. The mass-to-light ratio versus I-band central surface 
brightness for 12 galaxies from Bottema (1993) calculated accord- 
ing to equation (P3). Each point is marked by its NGC number 
(or G for the Galaxy). For the Galaxy Tj has not been divided by 
h. Squares are for inclined galaxies; triangles for face-on galaxies; 
circles are additional points for highly inclined galaxies where zq 
is measured (otherwise zg/R^ = 0.2 is assumed). The errorbars 
are estimated from the errors in the velocity dispersions. 



mass function (IMF) , particularly from the low-mass cut-off 
in the IMF. A few stellar population models are available 
(e.g., BerteUi et al 1994; Bruzual & Chariot 1993; Worthey 
1994). For a Salpeter IMF, the mass-to- light ratio derived 
by various authors appear to agree within an accuracy of 
25% (Chariot, Worthey & Bressan 1996). The predicted T/ 
depends on the metallicity and age of the stellar population. 
For a stellar population with age between 5-12 Gyr and with 
a solar metallicity, T/ is between 0.9-1.8 for a constant star 
formation rate (cf Table 3 in de Jong 1996). For an exponen- 
tial star formation law, the mass-to-light ratio is about 20% 
higher. The predicted values are in good agreement with the 
values derived from the instability analysis presented in this 
paper. In the comparison, we have neglected the uncertainty 
due to dust, since the TuUy-Fisher studies already attempt 
to correct for its effect. Furthermore de Jong (1996) argues 
that dust reddening probably plays a minor role in the color 
gradients in disk galaxies. Nevertheless, the dust correction 
remains a nuisance in these comparisons. 

A direct measurement of the mass-to-light ratio of ex- 
tragalactic disks requires detailed kinematic studies such as 
that described by Bottema (1993). Bottema (1997) derives a 
value of Ti = (1.7±0.5)/i after somewhat uncertain process- 
ing of his data. The measurement of T relies on the relation 
between vertical velocity dispersion a^, surface density E 
and vertical scale height zq. Essentially, the larger the value 
of E, the hotter a disk has to be at constant zq. For a disk 
which is well approximated by an isothermal sheet (cf. equa- 
tion |l|) 

o-^ = ttGEzo (21) 

(e.g. Binney & Tremaine 1987). Assuming an exponential 
radial profile we can write the mass-to-light ratio as 



T = 



2Gcr^Rd Rd 
Ld zo ' 



(22) 
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Figure 6. Shows em for the galaxies of Bottema (1993) given 
the values of T/ from Figure H 



Taking account of the possible variations in vertical density 
profile in the disk, equation (^2|) represents an upper limit on 
T (Wielen & Fuchs 1983). Bottema measures az[R = Q) for 
face-on galaxies, and asserts that aR{R — R^) = az{R~0) 
for inclined galaxies (this is true for the Galaxy). Only 4 
of the 12 galaxies are edge-on enough to measure zo, so an 
assumption also has to be made about its value in the other 
galaxies. Bottema discusses this problem in some depth. 
Where it can be measured directly zq/R^ ~ 0.2 and does not 
vary by more than about 20 percent (Barteldrees & Dettmar 
1994). 

Figure ^ shows T/ calculated from equation (^2|) (as- 
suming zo/Rd — 0.2 where no measurement is available) 
plotted against fio- The quantity plotted is independent of 
h. In the case of the Galaxy Tj has not been divided by 
h. The first impression of this figure is a marked trend of 
smaller T/ for higher surface brightness galaxies (low /io in 
magnitudes). The range of values of Tj is also large (about 
a factor of 10 between smallest and largest values). This is in 
marked conflict with the received wisdom regarding mass-to- 
light ratios of galactic disks, based on the uniformity of the 
colours of stellar populations (de Jong 1996). Indeed taken 
at face value, the trend in Figure ^ is so strong that when ap- 
plied to equation ( |l7| ) it should produce an anfi-correlation 
between e; and /io for constant md- 

How should we interpret Figure ^? We would like to see 
many more data before any strong conclusions are drawn. 
One can separate the galaxies by eye into two groups: one 
with high surface brightness (/io < 19.5) and low T/, and the 
other with lower surface brightness (/to > 19.5) and higher 
T/. Of the four galaxies with /to < 19.5, two have significant 
non-stellar contributions to their luminosity (3938 and 1566) 
which could bias them artificially towards the bottom left 
of Figure ^ Notice that few galaxies in the MF sample have 
surface brightness as high as these four galaxies. Taken on 
their own, the galaxies with /io > 19.5 have values of Tj 
which are perhaps not a strong function of /io. 

From equations (^) & we can obtain a simple ex- 
pression for Em in terms of kinematical quantities: 



Figure q shows the values of em for the galaxies in Figure 
^calculated from equation (p3|). They are all ^ 0.75 so the 
values of T/ in Figure ^ are consistent with the disks being 
globally stable. 

There remains a constraint on the average T/ as a func- 
tion of /io from the theory of disk formation. As was men- 
tioned in Section ^ the observed positive correlation in Fig- 
ure is consistent with A being at most a weak function of 
/io, particularly for the high surface brightness galaxies. This 
conclusion is expected to be only weakly dependent on the 
detailed model of disk formation — only the normalisation 
depends somewhat on the cosmological model (Mo, Mao & 
White 1997). The mass-to-light ratio T/ required to match 
the TuUy-Fisher relation is in the range of 1 to 2h, consistent 
with those obtained from the stability analysis presented 
here and models of stellar population synthesis. Whether 
it is actually consistent with Figure |^ is unclear until we 
have more data. If such a trend exists, the small scatter in 
the TuUy-Fisher amplitude A would require a nearly linear 
correlation between md and T/ (cf equation |l^) . 



7 DISCUSSION 

Global instability in a disk leads to bar formation, so it is 
natural to assume that barred galaxies form in this way (e.g. 
Sellwood 1996). If the correlation of e; with /io really refiects 
a correlation of em with /lo, then we should expect galaxies 
with smaller values of e; to form bars. No such effect is visible 
in the MF data. The distribution of barred types is indis- 
tinguishable from that of disk galaxies as a whole (Figure 
|l]) . Bar formation does not have to be spontaneous: it could 
be induced by a perturbation from a close encounter with 
another galaxy (Noguchi 1996), or from an interaction be- 
tween the disk and the halo. Perhaps all disks have T/ such 
that a bar does not form spontaneously, but they are equally 
susceptible to induced bar formation. This is consistent with 
the values of em in Figure ^ Note however that the barred 
fraction in the MF data is low (« 10%) compared to that 
reported by other authors (e.g. Sellwood & Wilkinson 1993 
give a fraction of around 30%). Further discussion of this 
important question is outside the scope of the present work. 

Clearly it is important to extend the range of data anal- 
ysed in Section This requires a good HI rotation curve for 
each galaxy, and high quality spectroscopy at least along 
the major axis of the galaxy. Ideally one would choose the 
brightest members of a large pre-defined sample (such as 
that of Mathewson & Ford 1996) and follow them up with 
a high spatial resolution spectrograph. Kinematic informa- 
tion in more than one dimension is also advantageous since 
it removes some of the uncertainties in the deprojection of 
the velocity ellipsoid. A number of two dimensional spec- 
trographs are due to come on line shortly, and these may 
be well suited to the problem. The biggest uncertainty in 
the determination of T/ will remain that associated with 
the value of zq. Efforts should therefore be made to anal- 
yse as many edge on galaxies as possible to try to improve 
on existing determinations of the distribution of zo /Rd (e.g. 
Barteldrees & Dettmar 1994). 
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APPENDIX 

Here we describe how values of R^, no and L were derived 
from the published quantities of Mathewson & Ford (1996). 
The published data list total magnitudes, and face-on cor- 
rected isophotal quantities: average surface brightness and 
isophotal diameter at /.i = 23.5mag/sec'^. Assuming an ex- 
ponential disk we have a surface brightness profile 

R_ 

Rd' 

and for a bulge to disk ratio / we have (outside the bulge) 



fi{R) = /io exp(— q), a ; 



(24) 



L{R) = id (1 + / - exp(-Q)(Q + 1)) . 



(25) 



From equation (^5|) the average surface brightness inside ra- 
dius R is 



fj. = = (1 + / - exp(-a)(a -f- 1)) , 



(26) 



where we have used Ld ~ 27r//o^d- Combining equations 
(M) and (Gfl) we obtain 



(27) 



p. ^ exp(Q)(l + /)^(Q-H) 

Given / and fi/fi we can solve equation numerically to 
find a. Then we have Rd — R/a, and fio ~ ^exp(a). 

In Table |l] the data labelled MF is obtained by setting 
f — (no bulge). Those labelled MFD have a crude bulge 
subtraction applied as follows. If the galaxy is classified Sbc 
or Sc a value of / is initially chosen randomly from a uniform 
distrubution in the range [0, .2]. If it is Sa or Sb / is similarly 
chosen in the range [0, .4]. This covers about 80% of the 
galaxies, and for the remainder we simply set f — 0. For too 
large a value of /, equation (^^ does not have a solution 
for positive real a. If / is initially too large, it is reduced 
by a factor of 3/4. This factor is applied repeatedly until a 
solution is found. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science DTJ5X style file. 
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